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' Abstract 

• The scalar two-loop self-energy master diagram is studied in the case of arbitrary 

O , masses. Analytical results in terms of Lauricella- and Appell-functions are presented 

for the imaginary part. By using the dispersion relation a one-dimensional integral 
representation is derived. This representation uses only elementary functions and is 
^ . thus well suited for a numerical calculation of the master diagram. 



1 Introduction 

The recent years have led to a considerable progress concerning the precision of measure- 
ments of some parameters of the standard model. Outstanding examples are the determi- 
nation of the the Z-mass and -width Also a precise determination of the top-mass can 
be expected in the next years. As a consequence, very precise theoretical evaluations of 
the theory are more and more required, including calculations of two-loop corrections. 

This paper focuses on the evaluation of the so-called scalar master diagram, which be- 
longs to the class of two-loop self-energy diagrams. Analytical solutions, involving general- 
ized hypergeometric functions, have been presented for all other scalar two-loop self-energy 
diagrams in the general mass case @-@]. But no comparable results have so far been found 
for the master diagram, whereas for special mass configurations or in asymptotic parameter 
regions there exist many solutions In this paper we are concerned with the general 

mass case. 
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After an introduction of our notations in sect. 0, we investigate in sect. |^ the imagi- 
nary part of the master diagram. Analytical results are presented for the two-particle cut 
contributions and three-particle cut contributions. 

Concerning the full calculation of the master integral, it must be noted that for 
two- loop integrals which contain a massive three-particle cut, the results in general can- 
not be expressed in terms of generalized polylogarithms with algebraic functions of the 
parameters as arguments. Thus, a promising approach to two-loop integrals uses integral 
representations which involve only elementary functions in the integrand and can therefore 
easily be evaluated numerically. For the master diagram a two-dimensional integral repre- 
sentation has become famous, which involves only simple logarithms and square roots. 
This approach has been generalized to other two-loop self-energy integrals in ref. [pH) . 

In sect. ^ we present a one-dimensional integral representation which only involves 
elementary functions. The standard algorithms for the numerical evaluation of one-di- 
mensional integrals show very good convergence properties. Therefore the representation 
presented in this paper allows a very fast numerical calculation of the master diagram. 

In sect. ^ we present some numerical comparisons of our results with those of pH 
demonstrating the agreement of the results obtained with both methods. 



2 Notation and definitions 




Figure 1: The master diagram T12345. 



The master integral (fig. |l|) is in the general mass case determined by six independent 
parameters, i.e. five internal masses and the external squared momentum p^. The name 
'master diagram' has been introduced by Broadhurst and is due to the fact that a cancel- 
lation of propagators in this diagram generates all other fundamental two-loop self-energy 
diagrams. The integral Ti2345(p^; m^, ml, 777.3, mg) finite in four dimensions. We use 



the convention 



1 r ... ... 1 



{kf — mi){{ki — py — ml) 
1 

m + k2Y - nil){kl -nil){{k2+pf - ml) 



where all masses have an infinitely small negative imaginary part. 



The integral T12345 has branch cuts along the positive real axis. The discontinuity across 
these cuts is related to the imaginary part, 



ATi2345(P^;"^i) = (Ti2345(P^; 



(2) 



where denotes the mass parameters. The integral can be calculated from the disconti- 
nuity with the dispersion relation 



Tl2345(P ' ) = ^ 7 



ATi2345(-5;mf 



so 



n ■ ) 



(3) 



where sq denotes the lowest branch point of the function. According to the Cutkosky rules 
[0, |l^, the discontinuity of the master diagram can be decomposed into four parts, each 
corresponding to a cut diagram (fig. 0), 
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where AT,^2345 5 AT;^^2345 ^''^^ ^"^uub^ AT;^^2345 denote the two- and three-particle disconti- 
nuities, cf. fig. 0. 






Figure 2: The contributions to the discontinuity of the master diagram. 



Each cut gives a contribution to the dispersion integral, which will be denoted accord- 
ingly. 
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For each contribution, the lower limit of the integration is determined by the threshold 



(2a) 



which is associated with the cut, i.e. s,^ 
(m2 + 7713 + rriiT-, sf^^ = {mi +771^ + 7n^y. 



(mi + m2)^ 4^''^ = (m4 + 7n^Y, 4' 



(3a) 



3 Calculation of the imaginary part 

3.1 Explicit results for the two-particle discontinuities 

The two-particle cut contributions to the discontinuity are given according to the Cutkosky 
rules as 



ATS5 = -\ ! d^hd^k22TiQ{kl)6{kl-7n\)2'KQ{-kl+p'^)6{{h-pf 



771^ 



{{ki + k2Y —7n1 — ie){k2 —7n\ — ie){{k2 + pY — mg — ie) ' 



(6) 



TT^ J 
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(fcf — mf + it){{ki — pY — 7712 + + k2)'^ — + ie) ' 



(7) 



In the evaluation of the expressions (^) and the discontinuity of the one-loop self-energy 
in four dimensions, 



ABoip^; ml ml) = ^ J d^k 27re(fc°)5(P - m?)27re(p° - ((p - kf 



mr 



. ^A(p2,mf,mi) 

= 27rz-5^ e(p^ - (mi +m2)^) , (8) 

occurs, where the Kallen function A is defined by X{x,y,z) = {x — y — z)"^ — Ayz. The 
expressions (H) and (0) yield 
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d% 27re(A;?)(5(A;2 _ ^2)27rG(/ - - p)^ - 



and equivalently, 



X (Co(p^ fc?, (fci - p)^; ml, ml, mlf^ 
= ABo{p^;ml,ml) (^Co{p^,ml,ml;ml,ml,ml)) , 

^^12345 = ^Bq{p^; ml, ml)Co{p'^, ml, ml; ml, m\, m^) , 



(9) 



(10) 



where Cq denotes the integral associated with the one- loop triangle diagram (fig. whose 




Figure 3: The triangle diagram. 



AT12345 cire hence given by a product of the one-loop self-energy discontinuity ABq and 
the one-loop triangle integral Co or its complex conjugate Cg. 

3.2 Integral representation for the three-particle discontinuity 

For the three-particle discontinuities AT12345 the Cutkosky rules yield 

Arn(3a) / 2 2 2 2 2 2\ 

ATI2AAV ■,m^,m2,m^,m^,mr,) 

= -—[ d% d% 27re(-A;° + /)5((A;i - p)^ - ml)27re{k^^ + A;0)5((A;i + ^2)' - ml) 

. 1 



analytical solution involves dilogarithms [0]. The two-particle discontinuities AT12345 ^-^d 



The contribution AT12345 requires no separate evaluation. It is given by 

/\rTi{3b) / 2 2 2 2 2 2\ A m(3a) / 2 2 2 2 2 2\ /io\ 

With the dispersion representation of the propagator, 
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1 f ^ 6{ki-t) 



dt ' \ , (13) 



/ci — mf + ie i t — m'f + ie 

— oo 

the two momentum integrations in (^) can be performed separately and yield 



.(3a) / 2. 2n _ 1 f , AEo(p^t,m|) ACo(t,mi,p2;m|,m|,m^; 



2m J t — mf + le 

(1713+1x14)^ 

The lower limit of the integration follows from the threshold of the Ai?o-function (H), while 
the upper limit guarantees that p"^ > (m2 + + 7714)^, i.e. that p"^ is larger then the three- 
particle threshold. In expression (0), the function ACq denotes the discontinuity of the 
one- loop integral Cq (fig. |^) with respect to the parameter pf. 

This discontinuity is the main ingredient in the dispersion representation of Co, which 
will be used in sect. |. The function Cq can be calculated as 
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1 f ACo{s,pl pf,mlm lml) 
Pi — ie 
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where the lower integration limit Sq = {mi + 7712)^ is the threshold belonging to pi, ACq 
is the discontinuity associated with the normal threshold and Coan is a contribution of an 
anomalous threshold. The discontinuity ACq is in general given by 

\ ri ( 22 2 2 2\ 

ACQ{s,p2,p^\mi,m2.,m^) 
= d% 27re(A:°)5(A;2 - m\) 27re(-A;° - p° - p°)5((A; + p2 + Pzf - ml) 



{k + pzY — ml — ie 



. log(a + h-ie)- \og{a - b - ie) ^ 
-2m — = Q{s - (mi + ms) ) (16) 

'{T.p)^ ~ s + ie J [Ap)^ — s + ie 



2tanh ^(b/a) ^, , ^9^ , 

-2m ) ' Q{s - (mi + m^)^) , (17) 

'(Jlp)^ — s + ie J [Ap)^ — s + ie 



where 



a{s,pl,pl;ml,ml,m^] 



= s{s + 2m\ -p\-p\-m\- mj) + {pj - pDimj - mj) , (18) 
b {s, pI,pI; ml, mj) = (T^p)^ - s + ie^ (Ap)^ - s + ie^J A(s, mf , ml) , (19) 

= \P2\ + \P3\ , Ap = \p2\ - \P2\ , 



and it is assumed that pi > and p| > 0. 

In order to produce a reliable numerical program for the master diagram, a careful 
discussion of the analytical properties of Cq is required. The Landau equations show that 
a leading singularity and thus an anomalous threshold occurs in Co if 

2 / 2 mipl + m2pl , 

TTin < nin.f.^ = mim2 . (20) 

■^""^ mi + 1712 

The branch point of this anomalous threshold is given by 

_ pl+pl + ml + ml-ml {pi - m^)(p| - ml) + I {pi ml ml) I {pi ml m|) 
^'"2 2 mi ^ 



with / (pj , mj,m,^) := ^pf — {mj — mfc)^ + i^yp'i — {f^j + "^fc)^ + • 

One can understand the effect of the anomalous threshold by an inspection of the behavior 
of the branch point Si as a function of ml taking into account that this parameter has an 
infinitely small imaginary part. In terms of the dispersion representation (|15D , in which 
Coan vanishes for m\ > mlf^^., the threshold si is a singularity of the logarithm in ACq as 
function of s. At ml = mlf^^ this singularity is equal to Sq. For m^ < mlf^^. the singularity 
Si crosses in some cases the integration contour, which consequently has to be deformed. 
The additional contribution Coan is therefore given by an integration over the discontinuity 



of the logarithm |T5 



Coan — 7^ / ;. , Z If ' (^2) 



1 r (27ri)2 ds 

which evaluates to square roots and logarithms 

ds 



I = = , ={ ^og{p\ - s + it) 

\JKs.pIpI){s-pI) \/x{pIpIpI) 

- log ( - s{pI +pI- pI) + {pI - pI)^ - p\{pI + pI) 

+ ^(Sp)2 - s + it^{Ap)^ - s + it^\{plplpl) + ^e)} . (23) 



From equation ( ^21) it can be observed that for si real and Si > Sq no separate treatment 
of C^an is required, but the discontinuity of the logarithm can be included as additional 
term to the logarithm in the discontinuity ACq (P^. 
Thus, one can distinguish the following six cases: 

• First case: If m\ > ml^^, no anomalous threshold occurs. 

• If ml < mli^j. the anomalous threshold has to be taken into account. The behavior 
of its branch point Si at m|^^^ depends on the value of X{pl ml mlf^^). 

— If X{pI ml mlf^j.) < 0, which is equivalent to \{pl ml mlf^^.) < 0, the anoma- 
lous threshold deforms the integration contour. Three subcases may occur: 

* Second case: Si is real and Si < Sq. 

* Third case: Si takes on complex values. 

* Fourth case: Si is real and Si > Sq. 



— If A(p2; '^i? ""^itftr) > ^he anomalous threshold deforms the integration contour 
only for complex valued Si. Two subcases can be distinguished: 

* Fifth case: Si is real and Si > Sq. Then the sheets of the logarithms in ACq 
change. 

* Sixth case: Si takes on complex values. 

For a numerical evaluation of the dispersion integral the following procedure for handling 
these cases can be applied. In the first case the discontinuity ACq can be calculated as 
given in (|l6l). In the other cases the sheets of the logarithms in ACq differ from ([16|) . In 
the second, third and sixth case the logarithms in (|1^) have to be replaced by 



log(a + b + ie) — log(a — b — ie) . 

In the fourth case the additional discontinuity of the logarithm yields 

J log(a + b + ie) — log(a — b — ie) — 2Txi for s < Si , 
1 log(a + b + ie) — log(a — b — ie) for s > Si . 

In the fifth case the correct expressions are 

log(a + b — ie) — log(a — b — ie) for s < Si and real b , 
log(a + b) — log(a — 6) — 2Tii for s < Si and imaginary b , 

log(a + b + ie) — log(a — b — ie) for s > Si . 

In the second, third and sixth case the contribution Coan of the anomalous threshold. 



given by the expressions (|22D and (|23D , has to be added. The argument of the second 
logarithm in (^3j) can get a negative imaginary part if Si has a nonvanishing imaginary 



part, i.e. in the third and sixth case. It has then to be replaced by the appropriate analytic 
continuation, which is given by 

- log ( - si{pl +pI- pI) + {pI - pI)^ - pUpI + pI) 



+ - + ze^iAp)^ -s,+ ie^\{plplpl)) - 2vr^ . (24) 

The expressions (^) and ( p4D are valid for p\ > {\p2 \ + IpsI)^- An analytic continuation in 
terms of pf is easy, but is not required in the context of this paper. 

We now return to the master diagram. Formula ([T^), which Broadhurst |^ took as 
starting point, is particularly useful, because the square root belonging to ABq cancels one 
of the square roots of ACq, which results in 
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AtSs (p^; m\) = — dt , , 2 tanh-^ . ' 2 I2 I2 I2 1,2^ > (25) 



p2 J t — m\ + ie a{t,m\,p'^,m\,m'l,m^) 

(7713+7714)2 



with a and b defined as in (|I^) and (|19|). The argument of the tanh is a rational function 

of t, p^, m\, m|, m\, m\, and of J \{t,p'^, m2)X{t, ml, ml). 



3.3 Representation of the three-particle discontinuities in terms 
of complete logarithmic elliptic integrals 

A partial integration in (^) yields for the three-particle discontinuity 
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x| (tanlr' -r'T'^ ) ■ (26) 

at \ a [t, 1712, p'^, ml, 171^,1714,) J 

Except for the logarithm, the form of ( p6D equals those of elliptic integrals. The zeros 



of the square root y A(t,p^, m2)A(t, mg, m^) as a function of t are located at (p±m2)^ and 
(m3 ± m^y. Therefore we introduce in the following the characteristic variables 

(1±± = {P± rn2Y - {rris ± 7714)^ . (27) 

In analogy to the case of the transformations of elliptic integrals to the Legendre normal 
form one substitutes 

(mg + m4)2g__ - (mg - m4)'^q_+x'^ 

which leads to 

^^12345 = 27ri / , (log(ao) + log(l - aix^) - log(l - a2X^)} 

{ ^(1-x2)(1-k2x2) ^ J 
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n {q— - cPg_+x2) 



2m I |log(ao) + log(l - aix^) - log(l - a2X^)| 

I ^ {I - X^){1 - K^X^) ^ ^ 



with 
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Ci C2 Cg C4 



rW + rWei" + 4"(c!")^ + rW(ei")'' + rW(4")^ ._ 

cf n(c<"-cf') tor,-l,...,4. 
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r4^'* = g g+_(m4 — ms) m3{p^ — — mi,) — 7714(1712 — — m 
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(^p2 _ ^2 _ ^2-j^^2 _ ^2 _ ^2^ _|_ 47723m4m5 ± A(m5, p^, ml)X{ml, m^, m^) 
Q^-* = (777,3 + "^4)Q'-+'?++(^2"^4 ^^3^4 — m^ml — m^ml — m^m^ + m^p^) , 

+2m4(p^ — ml)(—m2m3 + 2m^m2,m\ + ?)m\m\ 
^-m^m^p^ + 6m3?Ti4p^ + 5m4p^) 
+ (777-3 + ?T74)^(— 677727^73777,4 + 6/77277737774 — 27^7^7^74 + Am2m\ 

I 222 o42| 24 r)24 

+777277737774 — 277737774 + 77737774 — Zm^m^ 

— 277737774 + 7772^7^ — 277727773^7^ + T^gP^ 

+47/727/737774^^ + 2mlmlp'^ + m\m\p^ + Gm^m^p'^ 
—2m\p'^ — 2m\p'^ + 277737774/7'* — lOmlp"^ + P^)] 
+r775[ — 277737774(^7^ — 7772)^ 
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.4 I r,™2^2 I r,™2^2\ 
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rg^-* = Qm3m4{p'^ — ml + 7773 — 7774) (A(p^, 7774, 7775) — A(7772, 7773, 7775)) , 

fl^ Tr/*? \ / A A. 'y 9 

r\ = —2<. [7774 (p — 7772) (—277727773 + 77727774 + 4/77277737774 — 6/7727774 + 2//72//73P 

— 2/772//74P + 12/773/774^ — 10/r74P + /r74P ) 

+ (//73 — //74)^(— 6//72//73//74 + Gmlmlm^ — 2m\mi + lOmlp"^ — p^ 

A A 2 222io42 24 
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QO 0/1 0/1 
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/ \0//0 0\0 00 A 
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The integrals in (p9|) are either of the type of complete elliptic integrals or they have the 



structure 



LK(a, k) 
Ln(a, c, k) 



log(l — ax'^)dx 
(l-a;2)(l-K2a;2) 

log(l — ax^)dx 1 

(1-X2)(l-K2a;2) 1 - CX2 



(30) 
(31) 



The analogy of these functions to the elhptic integrals is obvious. According to our knowl- 
edge they have not been introduced in the mathematical literature. We call them complete 
logarithmic elliptic integrals. For sake of completeness one can introduce a third function, 



LE(a, k) 



Then all integrations of the type 



(1 — K^x'^) log(l — ax'^)dx 

^(l-x2)(l-/t2x2) 



(32) 



log(l — ax'^)dx 
I y^(l-a;2)(l-/«2a;2) 



r{x) 



(33) 



where r{x) is an arbitrary rational function, can be expressed in terms of LK, LE, LII, 
complete elliptic integrals, logarithms and polylogarithms. 

The functions LK, LE and LH are related to generalized hypergeometric functions of 
two and three variables, i.e. the Appell function Fi(q;; /?2; 7; a;, y) and the Lauricella 

function F2)^(a;/5i,^2,/?3;7;a;,l/) IT 
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The representation in terms of generalized hypergeometric functions helps to perform an- 
alytic continuations for LK, LK and LH. Series representations are given as 



LK(a, k) 
LE(a, K,) 
Lnfa, c, k) 



E 
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— (K^)™a" 
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(37) 
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The region of convergence is given by < 1, \a\ < 1, |c| < 1. 



In terms of these functions one can finally express the result for AT12345 as 
^TilUp'M) = 2vrz/.«| log(ao)(/^«K(«:) + ^/^«n(^c«,«:)) 

+ f 4'^LK(ai, K) + j2 hf^ Ln(ai, — cf \ 

V 1=1 1 — / 



(40) 



a result which involves in an analogous way complete elliptic integrals and the complete 
logarithmic integrals introduced in (|30|) , (|3T1) and (^). 

4 One-dimensional integral representations 

In many cases dispersion representations have proven to be useful for the calculation of 
two- loop self-energy integrals f^, ^. One can use this approach for further analytical 
calculations, or evaluate numerically the dispersion integral. In the case of the master 
diagram, the discontinuity in the general mass case is complicated to calculate, as shown 
in the last section. We will present two ways to simplify the evaluation of the dispersion 
integral. The first one is a realization of an idea of Broadhurst ||^. It simplifies the 
evaluation of the three-particle cut contributions by a partial integration in the dispersion 
integral. The second way relies on interchanging the integrations in both the two-particle 
cut contributions and the three-particle cut contributions of the master diagram and leads 
to a one- dimensional integral representation involving only elementary functions. 

4.1 Evaluation of the three-particle cut contribution 

Formula (|]) expresses the two-particle cut contributions to the discontinuity in terms of 
well-known functions, though involving dilogarithms. However, for the three-particle cuts, 
equation (p^ ) gives only an integral representation and expression ( P0| ) is a solution in 



terms of rather complicated series representations. 

As Broadhurst has pointed out |]^, it is useful to perform a partial integration in the 
dispersion integral for the three-particle cut contribution. One is then led to a suitable 
derivative, a'{p'^;m]), of the discontinuity. 

The function T can then be calculated from a' as 

1 °° 

T(p^; mf) = / ds log(s - - ie) a'{s] ml) . (42) 

p J 

so 

Equation (pSD yields for cr*^^")'(p^; mf ) of the three-particle cut contribution T^2U5 

a(-)'(pV?)=2 / -^#,ftanh-^ ^^'^ff^^-f^^^ - (43) 
J t — m{ dp^ \ a [t, 1712, p^, 711^,711^,711^) J 



After a lengthy decomposition into partial fractions one finally arrives at a surprisingly 
simple result, involving only three H-functions, i.e. complete elliptic integrals of the third 
kind, 

^ J 3 (2) 2 4 I (2) 2 2 

^il-x^){l-K^x^) ^=l ll{q-~-q-+cf^x^) 

i=l 

= / ; h^'^j: 

= /.(^)E/.rn(^cr,.), (44) 

with variables, 

2 _ q+-q-+ (2) _ ("^3 - m^f - mj 
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.(2) 
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q++q ' ^ (ms + m^y — mf ' 

(p^ — ml — m\){m2 — m\ — mf) — Am^m^rril ± yj\(ml,p'^, m1)X{ml, ml, m|) 
(p2 — ml — ml){m2 — m\ — m^) + Am^mj^ml ± \(ml,p'^, m1)\(ml, m\, m|) 

32m3m4 

y/q++q [^3(p^ — ""^4 ~ ^i) + ?Ti4(m2 — JTig — ml)]^ [(777,3 + ra/^Y — 7n,f] 



(2) , (2) (2) 

Tq^'' = (r7i2 — — 2m^m^{m\ — m^ + ml{ml + rrig) 

— 777.5(7772 + (m^ + m^Y) — ^7^(7772 + 7^73 — ml) , 
rj^"* = — rQ^''(7774 —7774) . 



4.2 An efficient one-dimensional integral representation for the 
master diagram 

The dispersion integral for the master diagram involves in the integrand dilogarithms for 
the two-particle cut contributions and complete elliptic integrals for the three-particle cut 
contributions (|1^), as presented in the previous section. In this section we use a different 
approach in order to derive a one- dimensional integral representation with elementary 
functions. This is achieved in the following way. 

For an evaluation of the two-particle cut contributions it proves useful to introduce the 
dispersion representation of the function Co as given in equation (|T^). After an insertion 
of (^ and ([T5| ) into (^, the integrations in the dispersion integral can be interchanged and 
lead to 

rr(2a) / 2 2 2 2 2 2\ 
Tl2Ub\P ;"^1,"^2'"^3'"^4)"^5) 

so 



2niJ J J (s — p"^ — ie)(t — s + ie) 
to «o 

1 7 ABo{s;mj,ml) , 2222 2^y 

+— / as : [Coan{s,m^,m2;m^,m^,m2)} 

2m J s — — le ^ ' 



so 

00 



1 1 ^^ (ACo(t.m;.ml;ml.mimi))- ^ _ ^ 

2m J t — p^ — le ^ ^ 

(m4+m5)2 



,1 1 A ^Bo{s-ml,ml) ( 2222 2^^* u<\ 

+— / rfs 5 ^ Coan(s,mi,m2;m5,m4,m3) . (45) 

2m J s — — 26 ^ ^ 

(r)ii+m2)2 

Thus, one ends up with two one-dimensional integrals. The infinite and the constant parts 
of the two i?o-functions |n| in (^5]) cancel. In both integrals the integrands are composed 
of logarithms and square roots only, as can be seen from equations ([T6|) and (^31) . The two 
integrals can therefore easily be evaluated numerically. 
For AT]-2345 the analogous calculation yields 



rr{2h)f 2 2 2 2 2 

1^ '[p ■,mi,m2,m^,m^,m 



^ 7 ^ACo{t,mlml;ml,mj,ml) /„ . 2 2 2x /, 2 2 

— / dt ; [Bo{p ; m^) - [Bo{t; m^, 

2m J t — p'^ — le ^ ^ 

(mi+m2)^ 



1 7 /\Bo{s]m%ml) 2 2 2 2 2n 

+ — / as : Coan[s,m^,m^]m2,m^,m^) . (46) 

2m J s — p^ — le 

The three-particle discontinuity (|Tip can be transformed with the partial integration 
6]), yielding 



ATil^Us;m^) = 2m J dt log {-^^ - 1 + ze 



{m3+m4)2 



J A(s, t, ml) R{s, t, ml, mj, ml, ml) 

X ^ , (47) 

s 

where R{s) is a rational function of s. In the general mass case, i.e. if all masses are 
different, it has only first order poles. This form is well suited for evaluating the dispersion 
integral, 

~ 2m 

so 



^12345 P ^^i ) = t;—- ds : , 4^ 

2m J s — — ze 

so 

because R{s) can be decomposed into partial fractions. 



^\it,mlml) 1 ( rf^g + rf ^ (3) (3)' 

tmi {s - Si){s - S2) \ (s - Ss){s - S4) 



^ Xjt, ml ml) ^ f^{t, ml m\, m\, m\) ^^^^ 



tml ^ s- Si 



(3) 

The values of the and Si will be presented below. An interchange of integrations, 



oo (v^-m2)^ oo oo 

ds J dt= J dt J ds, (50) 

(m.2+m.3+m4)2 {ms+rm)^ (7713+014)2 (ma+v^)^ 



/ 



leads to integrations of the type 

ds \J\{s,t,ml) 

s 

(Vi+m2)2 



/ ^ = B,{s,,t,ml) , (51) 

J S — Si s 



i.e. one-loop self-energy integrals. In the final expression all ultraviolet divergencies of 
these -Bo-functions cancel. The result for the three-particle cut contribution of the master 
diagram is then 



T^^^"^) (p2; m^i) = J dt log! + 



X 



\{t, ml mi) 4 Q^^p2. ^2) _ 5^(3.. 



'y^ f - ^ 1 "^2/ -^uy-^ii '"-2; (52) 



7 = 1 

with 



tml ^ - p2 _ 



Sl/2 



t + ml + ml + ml~Tnl {t - ml) {ml - ml) ± ^A(t, m§, ml)A(mi, m^, m^) 



2 2mi 



S3/4 = ("^2 ± Vi) 



2 



/i = ^ + , /2 = /i(si ^ S2) 

n (^1 - ^7) 

i=2 

/3 = -j=j-7 r, J4 = JslSs ^ S4) , 

n (S3 - Si) 

77^3 

(3) 1/0 2 2 I 2\ 2/2 2n 

r| = t[2m2 — m^ + m^) — m2[m^ — m^) , 

= - "^2)(^("^5 - ^3) - ml{ml - mj)) , 
^,^_ml{ml-t-ml) 



X{t,mlml) 

(3) 2 

t{2m\m\ — m\m\ + m\m\ + m\m\) + m\{m1 — m\){m\ — m|) + m\m\{m\ — m|) 



+ 



A(t,ml,mi) 



In order to obtain the second three-particle cut contribution T^^^\ the masses have to 
be interchanged according to mi m2 and m4 ^ ms. 

In some parameter regions one has to be careful when evaluating the i?o-functions in 
(|52D. The external momentum variable can be assumed to have an infinitely small 



positive imaginary part. Consequently, BqIjP' + ie;t,m\) can be evaluated in the usual 
way. The functions i?o(s3/4; t, m^) are always evaluated below the threshold, and present 



therefore no difficulty. However, care has to be taken in the evaluation of 50(51/2; t, m^). A 
brief inspection of the asymptotic behavior of Si/2 shows that the functions Bq{si/2; t, m^) 
have to be evaluated below the threshold for large values of m\. In fact, no difficulty occurs 
for m\ > (m2 + ms)^. But Si and S2 get an imaginary part for 

(m2 - ms)^ < ml < {1112 + m^Y . 

In that case, the functions Bq{si/2', ml, m^) have to be evaluated for complex values of the 
variables Si/2, which is not too complicated. For m^ < {m2 — fn^Y, the imaginary part of 
the variables si/2 vanishes again. Then one can in general evaluate the i?o-functions with 
the usual algorithm which assumes -Bo(si/2 + ie^t^m^)- However, there exists a threshold 



tthri = j-^(x{ml,ml,ml) +2ml{ml + ml) 

+ {ml - ml - ml)^ \{ml,ml,ml) + AmlmU , (53) 



and for t > tthri the function Bo{s2;t,ml) has to be evaluated as Bo{s2 — ie^t^ml) = 
{Bo{s2 + ie;t,ml)y. 

Another complication occurs in that case, if S2{t) equals at a value t > tthri- Then 
one denominator in (|52D vanishes at a point where the numerator takes on the value of the 
discontinuity of a i?o-function. This happens at t = tthr2, if the following conditions apply 

X{p^mlml)>0, (54) 
tthri < tthr2 = (ml{p'^ + ml + ml + ml- ml) + {mj - ml){ml - p^) 

+ \/ ^("^2, ml)\(y, ml, ml)^ , (55) 
tthr2 >to = (mg + m4f , (56) 
mj + ml-mj- \J \{ml,ml,ml) > . (57) 



The condition (Q) is necessary for tthr2 to be real valued. The conditions (|55| ) and ( [56| ) are 
obviously necessary for the complication in question. The condition (^) can be understood 
from the asymptotic behavior of S2 for large values of t. Thus, if m| < (m2 — m^^ and 
the conditions ( ^^ - |57| ) apply, the principal value and the residue of the integral have to be 
calculated. 

The calculation of the residue involves the derivative of S2 with respect to t. 



dso 1 



2 2 I 2 I / 2 I 2 ±\ 

3 — m2 + mg + (mg + m^ — t)> 



— - - — ^ , m, - m2 + + [m^ + m^-t . — — 5 ^ 

dt 2ml \ ^ b \ i 4 X{t,ml,ml) 

The calculation then yields 



(5^ 



J S2 - 



to 



Some special mass cases have to be considered. Double poles occur in R{s) if m2 = 0, 
in which case S3 = S4, or if X{m2,ml,'m'^) = 0, in which case Si = S2- In the case of 
ml = the partial integration (p6D is modified. For T^^''\ these conditions amount to 
ml = 0, X{ml, m\, m|) = or m^ = 0. The decomposition into partial fractions leads then 
to modified results involving functions 



ds y^X{s,t,m^) dBois,;t,m^) 

- • (62) 



(s — Si — ie^ s dsi 

(v^+m)2 

These derivatives of the i?o-function are finite in four dimensions. Another modification 
and considerable simplification of the decomposition into partial fractions occurs if m\ = 0. 

In our evaluation of the three-particle cut contributions the symmetry of the master 
diagram with respect to an simultaneous interchange of the masses ml m| and ^ m^f 
is hidden. Broadhurst has evaluated the master diagram for many cases of physical interest 
0. Note that all these cases can be understood to belong to the special cases specified 
above, if one takes the symmetries of the master diagram into account. 



5 Numerical comparisons 

The formulae (^5|) and (^2]) provide an efficient algorithm for a numerical calculation of 



the master diagram in the general mass case. The results have been checked by a compar- 
ison with Kreimer's two-dimensional integral representation [0]. We implemented both 
algorithms with C++ and performed a test for more than 150000 different values of the 
parameters. This should be sufficient to check the parameter region. 

For the one-dimensional integral representation an adaptive Gauss-Kronrod algorithm, 
implemented in the QUADPACK-routines |T3] proved suitable. For the two-dimensional 



integral we did not succeed with a Gauss integration. But following ref. |]TT1, we used the 



VEGAS-routine |]T9[ which implements an adaptive Monte Carlo integration algorithm. 

The following tables are intended to compare typical results of the two algorithms. 
We have chosen a parameter region where has a similar magnitude as the internal 
masses, where the integral representation is most useful to handle. Table |I| gives the results 
of the one- dimensional integral representation, table |^ of the two-dimensional integral 
representation. In each row the results for the real and the imaginary part, the reported 
error and the CPU-time are presented. To get a better estimate of the real errors and some 
impression about the convergence qualities of the algorithms, two runs have been performed 
for each parameter set. In case of the one-dimensional representation the required accuracy 
has been 10~^ for the first run and 5 x 10~^ for the second run. In the case of the two- 
dimensional representations we used 10^ function evaluations for the first run and 5 x 10^ 
function evaluations for the second run. The numerical results for the imaginary parts are 



reported even in those regions, where is below all thresholds, because the magnitude of 
the imaginary part gives another estimate for the errors. 

The precision of the one-dimensional integral representation can be increased further by 
introducing asymptotic formulae for the integrand at large values of the integration vari- 
able. With one-dimensional integral representations for other two-loop self-energy integrals 

^ we made the experience that simple asymptotic formulae increase the maximal preci- 
sion by several orders of magnitude. We are going to implement this also for the algorithm 
presented in this paper. A first version of the program can be requested from the authors. 

The programs have been run on a workstation DEC 3000 AXP. Note that the error 
reported by the VEGAS-routine gives the size of a standard deviation, while the error 
reported by the QUADPACK-routine is not so clear defined, but is in general highly over- 
estimated. This can be observed from the tables by a comparison of the results with 
different precision. 
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B 


C 


D 


E 


0.1 


< 


-2.87238e - 


01 


-1.84514e- 


09 


4.9e- 


05 


0.6 


0.1 


< 


-2.87238e - 


01 


-1.57367e- 


10 


2.5e- 


06 


0.8 


0.5 


< 


-2.94592e - 


01 


-1.89266e- 


09 


5.0e- 


05 


0.7 


0.5 


< 


-2.94592e - 


01 


-1.61303e- 


10 


2.5e- 


06 


0.8 


1.0 


< 


-3.04522e - 


01 


-1.95521e- 


09 


5.0e- 


05 


0.7 


1.0 


< 


-3.04521e - 


01 


-1.66484e- 


10 


2.5e- 


06 


0.8 


5.0 
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-4.52521e - 


01 


-2.62945e - 


09 


5.6e — 


05 


0.7 


5.0 


< 


-4.52520e - 


01 


-2.22337e - 


10 


2.8e- 


06 


0.8 


10.0 


> 


-4.88154e - 


01 


-3.53218e- 


01 


4.9e- 


05 


0.7 


10.0 


> 


-4.88153e - 


01 


-3.53217e- 


01 


2.3e- 


06 


0.9 


50.0 


> 


1.73902e- 


01 


-1.18080e- 


01 


4.3e- 


05 


1.4 


50.0 


> 


1.73901e- 


01 


-1.18080e- 


01 


2.2e- 


06 


1.8 



Table 1: The master integral Ti2345(p^; mf, m|, 77^4, rrig), calculated with the one- 
dimensional integral representation. The masses are m\ = 1, = 2, m| = 3, ml = 4, 
mf = 5. A: '<' indicates, that is below all thresholds, i.e. the imaginary part is zero. B: 
Numerical result for the real part. C: Numerical result for the imaginary part. D: Error 
reported by the QUADPACK-routine. E: CPU-time in seconds. 

As expected, the one-dimensional integral representation is considerably faster. The 
difference is more important if higher accuracy is required. Note that the CPU-time for 
the one-dimensional integration is nearly the same for the two runs for all values of p"^, 
while the accuracy has increased by a factor of 20. For the two-dimensional integration, 
on the other hand, a difference in the computing-time by a factor 5 causes about the same 
increase in accuracy. While the one-dimensional integration gives results accurate to about 
6 digits (regarding the real error) in about 1 second, the two-dimensional integration needs 
about 85 seconds to achieve an accuracy of between 3 and 4 digits, depending on the 
parameters. 
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-2 94846e - 


01 


4 17190e - 
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1.6e - 
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01 


-4.88930e - 
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04 


8.3e- 
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01 


-9.65061e - 


05 
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10.0 
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-4.88270e - 


01 


-3.52664e - 


01 


4.8e- 
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-4.88250e - 


01 


-3.53326e - 


01 
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04 


86 


50.0 


> 


1.73864e- 


01 


-1.18080e- 


01 


6.9e- 


05 


17 


50.0 
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1.73907e - 


01 


-1.18072e- 


01 


1.5e- 


05 


85 



Table 2: The master integral Ti2345(p^; mj, 7713, m|, m^), calculated with the two- 
dimensional integral representation. The masses are vn\ = 1, ml = 2, ml = 3, ml = 4, 
= 5. A: '<' indicates, that p"^ is below all thresholds, i.e. the imaginary part is zero. B: 
Numerical result for the real part. C: Numerical result for the imaginary part. D: Error 
reported by the VEGAS-routine. E: CPU-time in seconds. 



6 Conclusion 

For scalar two-loop self-energy integrals analytic expressions have been presented in the 
literature for many special cases P, ^. Asymptotic expansions and Taylor series are another 
important approach to these integrals, see e.g. [§. Furthermore, for diagrams containing a 
self-energy subloop, series representations have been presented which are analytical 

results for the general mass case. No comparable analytical result has so far been found for 
the master diagram in the general mass case. Anyhow, series representations always have a 
restricted region of convergence, i.e. in each particular region up to the closest threshold or 
pseudo-threshold. In the general case of the master diagram with its many parameters this 
will lead to considerable complications concerning representations in terms of generalized 
hypergeometric functions or in terms of Taylor series. 

In this paper we presented analytical results for the imaginary part of the master 
diagram T12345 and a one-dimensional integral representation for T12345. In all our results 
the two- and three-particle cut contributions are separated, which gives some additional 
insight into the analytical structure of the function T12345. Together with corresponding 
results P, ffl for the other fundamental two-loop self-energy integrals and with a tensor 



reduction formalism for the two-loop self-energy integrals ||2^, efficient algorithms are 
now available for the calculation of two-loop self-energies in the general mass case. The 
numerical implementations should be sufficient to perform complete calculations in the 
standard model at the two-loop level, e.g. to evaluate the fermion or gauge-boson self- 
energies. 
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